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ARTICLE INFO ABSTRACT 
Article History: BACKGROUND AND OBJECTIVES: The estimate of reserved claims is based on the 
Received 24 July 2023 prediction of the final amount of claims that have not yet been settled and which 
Revised 08 November 2023 the insurer has undertaken to pay. The micro-level approach estimates the amount of 
Accepted 18 December 2023 each unsettled claim separately. In this context, the time taken to settle a claim is an 


important variable, as large claims usually take longer to settle and may not be paid 
in the relevant financial period. Non-payment of claims may be due to a lack of timely 


Keywords: reporting, a high volume of cases or legal complications. Therefore, there may also 
Censored data be censored data. In this article, we use the copula function approach to model the 
Copula function dependency structure of the settlement duration and claim amount variables. 
Dependence structure 


METHODS: In this study, the amount of each unpaid claim is estimated separately using 
a micro-level approach. For this purpose, the duration of the settlement of each claim 
is considered as a variable that depends on the amount of the claim. Based on the 
modeling of the dependency structure between the claim amount and the settlement 
duration using the copula function and using the general characteristics of the claims 
as predictor variables, each unpaid claim is estimated. The marginal distributions of 
the claim amount and the settlement period based on covariates explain the stochastic 
behavior of each of them and are modeled separately. In the copula function approach, 
the dependence structure of the variables can be considered separately from their 
marginal distributions. By choosing or constructing an acceptable and appropriate 
copula function with the dependence structure of the claim amount and the duration 
of its settlement, and taking into account the particular characteristics and conditions 
of inflation, the amount of an unpaid claim can be estimated more accurately. In this 
study, the accuracy of the estimation and the validity of the proposed model are 
evaluated using simulations. 
FINDINGS: The proposed method is implemented for the data collection of one of the 
insurance companies regarding the cases of the employer’s professional liability claims 
in 8 years, from March 2013 to March 2021, which includes unpaid claims. In order to 
*Corresponding Author: evaluate the error in estimating the claim amount with this method in this data set, 
1000 samples are taken from the data whose claim amount is known and settled and 
each time a number similar to the actual number is censored, the known values are 
censored and their amount is estimated. Then the estimation error is calculated using 
ORCID: 0000-0002-8831-3372 the criteria mentioned in the article. And finally, it can be seen that the results have a 
good accuracy rate. 
CONCLUSION: In the article, it can be seen that the selected copula function has 
DOI: 10.22056/ijir.2024.03.06 acceptable results for the estimation of reserved claims. 
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پیشینه و اهداف: برآورد ذخيرة خسارت‌های معوق براساس پیش‌بینی مبلغ نهايى خسارتهايى است که هنوز تسويه 
نشده است و بیمه‌گر متعهد به يرداخت آن است. در رويكرد سطح خرد. ميزان هر خسارت تسويهنشده به‌طور جداگانه 
برآورد می‌شود. در اين راستاء مدت‌زمان تسوية هر خسارت متغیر مهمی است. زیرا تسوية خسارات بزرگ معمولا بیشتر 
طول می کشد و ممکن است در دورة مالی مربوطه پرداخت نشود. تسویه نشدن مطالبات ممکن است به‌دلیل گزارش نشدن 
به‌موقع. حجم CYL‏ پرونده‌ها و يا عوارض قانونی باشد. ازاينروء ممکن است داده‌های سانسورشده نیز وجود داشته باشند. 
ما در اين مقاله از روبکرد تابع مفصل برای مدل‌سازی ساختار وابستگی متغیرهای مدت‌زمان پرداخت خسارت و مبلغ 
خسارت. استفاده می‌کنیم. 

روش‌شناسی: در اين تحقیق, با رویکرد سطح خرد. میزان هر خسارت تسویه‌نشده به‌طور جداگانه برآورد می‌شود. برای اين 
منظورء مدت‌زمان تسوية هر خسارت به‌عنوان متغیر وابسته به میزان خسارت در نظر گرفته می‌شود. براساس مدل‌سازی 
ساختار وابستگی بين مبلغ خسارت و مدت تسويه با استفاده از تابح مفصل» و همچنین با استفاده از ویژگی‌های کلی 
خسارت‌ها به‌عنوان متغیرهای پیش‌بین ی کننده. هر خسارت تسویه‌نشده برآورد می‌شود. مدل‌سازی جداگانة توزیع‌های 
حاشیه‌ای مبلغ خسارت و مدت تسویه براساس متغیرهای کمکی. رفتار تصادفی هریک از آنها را توضیح می‌دهد. در رویکرد 
alli‏ مفصل ساختار وابستگی متغيرها را می‌توان جدا از توزیع‌های حاشیه‌ای آنها در نظر گرفت. با انتخاب يا ساخت تابع 
مقصل قابل قبول واا ماخر وایستگی مبلغ خسارت و منشزمان Bed gag lipid‏ كرفي ویزگی‌های ناش و 
شرایط تورم. می‌توان میزان هرخسارت پرداخت‌نشده‌ای را با دقت بیشتری تخمین زد. در اين تحقیق از شبیه‌سازی برای 
ارزیابی صحت برآورد و اعتبار مدل پیشنهادی استفاده می‌شود. 
بافته‌ها: روش پیشنهادی برای مجموعه داده‌های یکی از شرکت‌های بیمه مربوط به موارد اعلامی ناشی از مسئولیت 
حرفه‌ای کارفرما در A‏ سال که شامل مطالبات پرداخت‌نشده نیز هستند اجرا می‌شود. برای ارزیابی خطای برآورد میزان 
ذخيرة معوق با این روش در اين مجموعه داده. از بين داده‌هایی که مبلغ خسارت آنها معلوم است و تسویه شده‌اند به تعداد 
۰ بار نمونه‌گیری شده و هر بار تعدادی مشابه تعداد واقعى سانسورشده. مقادیر معلوم سانسور می‌شوند و ذخیرة معوق 
آنها برآورد می‌شود. سپس خطای برآورد با معیارهایی که در اين مقاله آمده محاسبه شده است. در نهایت مشاهده می‌شود 


که یافته‌ها ازدقت خوبی برخوردارند. 
نتبجه گیری: در مقاله ديده می‌شود که تابع rade‏ انتخاب‌شده نتایج قابل قبولی برای برآورد ذخيرة خسارت‌های معوق 


داشته است. 


توجه: مدت زمان بحث و انتقاد برای اين مقاله تا | اکتبر ۲۰۳۴ در وب‌سایت WIR‏ در «نمایش مقاله» باز است. 
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صديقه شمس و همکاران 


مقدمه 

موضوع SES‏ و محاسبة آن یکی از مسائل بسیار مهم شرکت‌های 
بیمه است و صنعت dar‏ کشور به‌دنبال روشی مناسب با شرایط 
بومی کشور برای حل اين مسئله است. تاکنون اين برآورد در سطح 
کلان و با روش نردبان زنجیری انجام شده است و ایرادهای زیادی 
ازسوی don‏ مرکزی به آن وارد بوده است. هدف اين مقاله ترویج و 
بومی‌سازی روش برآورد ذخيرة معوق در سطح خرد با استفاده از تابع 
مفصل است. 

برآورد ذخيرة خسارت معوق متکی به پیش‌بینی میزان نهایی 
خسارت‌هایی Cul‏ که هنوز تسویه نشده‌اند و بیمه‌گر متعهد به 
يرداخت آن است )2019 (Karimi,‏ ممکن است خسارت‌ها در 
دورة مالى مربوطه پرداخت نشونده زيرا معمولاً خسارت‌های بزرك: 
مدتزمان بيشترى براى تسويه شدن لازم دارند بنابراين با داده‌های 
سانسورشده مواجهيم (به‌عنوان يك منبع براى داده‌های سانسورشده 
می‌توان به )2003( Klein and Moeschberger‏ مراجعه كرد). در 
واقع» در تاریخی معین. يا خسارتی تسویه شده است که در اين 
صورت مقدار نهایی آن مشخص است. يا هنوز تسویه نشده است و در 
اين صورت ale‏ نهایی خسارت سانسور شده است. 

دو رویکرد کلی برای برآورد ذخيرة خسارت معوق وجود 
دارد: رویکرد سطح OU‏ و رویکرد سطح خرد. در رویکرد سطح 
کلان معمولا از روش نردبان زنجیری استفاده می‌شود که در 
Mack (1993) .Mack and Venter (2000) ۸۱۷/۵۲۲۵۱۱ (2000).‏ 
England and Verrall (2002)‏ و )2014( Saluz et al.‏ شرح داده 
شده است aS)‏ موضوع بحث اين پژوهش نیست). در اين مقاله از رویکرد 
سطح خرد استفاده می‌شود و برای پیش‌بینی مقدار خسارت معوق از 
مشخصات شناخته‌شدة خسارت به‌عنوان متغیرهای پیش‌بینی کننده 
استفاده می‌شود. برآورد در سطح خرد. در عمل کارایی خوبی دارد. اما 
از نظر اجرایی پیچیده است )1994 (Hesselager,‏ در رویکرد سطح 
خرد مبلغ هر خسارت تسویه‌نشده به‌صورت جداگانه برآورد می‌شود 
(Norberg, 1993; Norberg, 1999; Jin, 2017; Pigeon et al., 2014)‏ 
و مدت‌زمان پرداخت هر خسارت متغیر مهمی به شمار می‌رود 
(Yanez and Pigeon, 2021; Taha et al., 2021; Lopez‏ 
.and Milhaud, 2021)‏ در پژوهش‌هایی ath‏ رویکرد 
Ayuso and Santolino (2008), Jin (2017)‏ متغیرهای مدت زمان 
پرداخت خسارت و ale‏ خسارت به‌طور جداگانه بررسی شده‌اند و 
Antonio et al. (2016)‏ مدلی چندحالتی برای ale‏ خسارت در 

رویکرد al‏ مفصل برای تبیین ساختار وابستگی بين 
متغیرهای مدت‌زمان پرداخت خسارت و زمان وقوع خسارت را 
Zhao and Zhou )2010(‏ مطرح کرد. ما در اين مقاله از رویکرد 
Lopez et al. (2016)‏ برای مدل‌سازی ساختار وابستگی متغیرهای 
مدت‌زمان پرداخت خسارت و ale‏ خسارت استفاده م ىكنيم. در 
اين رویکرد اطلاعات در سطح خرد تحلیل می‌شوند. ابتدا مدل‌سازی 
(جداگانه) توزیع‌های حاشیه‌ای میزان خسارت و مدت‌زمان تسویه 


براساس متغیرهای کمکی انجام می‌شود. سپس اين توزیع‌های 
حاشیه‌ای با یک fade al‏ به یکدیگر پیوند داده می‌شوند و ذخيرة 
خسارت‌های معوق با استفاده از تابع مفصل برآورد مى شود. يس از 
انتخاب توزیع‌های حاشیه‌ای و توزیع توام تبیین می‌شود. در بخش 
سوم با استفاده از شبیه‌سازی کارایی مدل ساخته‌شده در بخش دوم 
ارزیابی می‌شود. در بخش چهارم مدل ساخته‌شده برای پیش‌بینی 
ذخيرة خسارت‌های معوق یکی از شرکت‌های بیمه به کار برده خواهد 


مبانی نظری پژوهش 

فرض كنيد متغیرهای تصادفی M‏ و 7 به‌ترتیب نشان‌دهندة 
مقدار خسارت و مدت‌زمان تسوية خسارت. از زمان اعلام گزارش 
باشند. همچنین فرض كنيد متغیر تصادفی D‏ زمان بين اعلام گزارش 
تا تاريخ شروع مطالعه. و بردار K‏ بردار متغیرهای بيش بين ی کننده 
باشند. در اين صورت تعریف می کنیم: 


Y =inf (T,D) 
۵ - وما‎ 
N=6M 


و برای تعداد 77 خسارت. مجموعه‌ای مستقل و هم توزیع 
Xi Nii ce,‏ رز رزة (Yi‏ از بردار تصادفی (Y, 5, D, ×, N)‏ به‌ازای 
i=l,...,n‏ هستند. در مدل سطح خرد. M,‏ مقدار خسارت 
شمارة 7 است. بردار p‏ بعدی» X, ER?‏ بردار متغیرهای کمکی 
كه خصوصیات هر خسارت است. T,‏ مدت‌زمان تسوية خسارت است 
(يعنى زمان بين اعلام وقوع خسارت و تسوية نهایی Col‏ اين متغیر 
ممكن است با D,‏ سانسور شود متغير ,8 .يك متغير نشائكر اسث 
كه برای موارد تسویه‌شده Gia]‏ و برای موارد تسويهنشده 0-:8 است. 
همچنین متغیرهای Y; =inf(T;, DJ‏ و6۸ = Ni‏ را تشكيل 
می‌دهيم. اين مدل مشابه مدلی است که در )2016( Lopez et al.‏ 
استفاده شده اسث» با اين تفاوت كه در اينجا فرض بر اين است كه 
D;‏ مشاهده شده است. 

فرض مىكنيم (T,M, X)‏ مستقل از D‏ باشند. به‌عبارت‌دیگر 
طول بازة زمانی اعلام خسارت تا شروع مطالعه. مستقل از متغیرهای 
کمکی X‏ مبلغ خسارت. طول دورة زمانی پرداخت خسارت. است. 

در بخش بعدی روش‌های برآورد توزیع‌های حاشیه‌ای 
- ۸ | 1 و M|X=x‏ شرح داده می‌شوند. 


برآورد توزیع‌های حاشیه‌ای 
برای مدل‌سازی توزیع شرطی × = × | ٠1‏ از تابع نرخ خطر 
dP(T > [t, t+ dt]|X = x)‏ 


h(t /x) = lim 0 


YYY 


ذخیره خسارت‌های معوق در سطح خرد 


و مدل رگرسیون Cox and Oakes (1984) (SIS‏ به‌صورت 
A(x) = ho) exp(alx)‏ که در hool‏ یک تابع ياية ناشناخته 
و a‏ بردار پارامتر نامعلوم asl‏ برآورد a‏ را می‌توان با استفاده از 
روش حداکثر شبه‌درست‌نمایی انجام داد» (در نظرية آماری تابع 
از هاش حا اسع كاريرت igs‏ 
تقریبی برای تابع احتمال مجموعه‌ای از مشاهده‌شده ارائه 
کند که ممکن Cal‏ روش محاسباتی ساده‌تری برای تخمین ارائه 
پارامترهای مدل ارائه کند). 

تابع خطر تجمعی ds‏ (5) رط Ho (= fy‏ را می‌توان با 

استفاده از برآوردگر پرسلو )1994( Burr‏ به‌صورت ناپارامتری برآورد 
کرد. 

مدل سازى توزیع شرطی M|X=x‏ براساس مدل خطی تعمیم 
یافته است )2010 ,۲۳۷۲6۵00 200 1972120560 «NelderandBaker,‏ 
در اين مدل یک y(M) bas‏ شناخته شده از M‏ در نظر 
i . yes‏ ۲ 7 
گرفته می‌شود و فرض می‌شود × (M ( /×[( = Bo‏ 9)۴]۷ 
كه در آن © يك تابع پیوند شناخته شده است. برآورد پارامتر By‏ 
را می‌توان با استفاده از 


= 


B = arg max, Zù, Wa toe fp (Y(ND|B Xi) 5 


T 7‏ ۱ 
كه در آن f88 X)‏ چگالی شرطی از y(M)|X‏ با پارامتر 
cul B‏ وزن‌ها 


0 
i nc) 


هستند که در آنها 9 = ۳ تابع بقا تجربی و یک 
براوردگر از S (T)=P(C2t)‏ است. رابطة O)‏ و (V)‏ مربوط به 
برآورد حداکثر درست‌نمایی وزنی است که وزن‌ها برای اصلاح کاهش 
اریبی ناشی از سانسور اعمال شده‌اند. 


مدل‌سازی براساس تابع مفصل 

در مدل سطح خرد. متغیر M pag T‏ وابسته هستند. در 
تسویه شدن. با مقدار بزرگی از M‏ همراه باشند. زیرا به احتمال زياد 
مطالبات کوچک به‌سرعت تسویه می‌شوند. مزیت مدل‌سازی ساختار 
وابستگی متغیر 1 و متغیر M‏ از دیدگاه تابح مفصل در اين است که 
می‌توان ساختار وابستگی را به‌صورت جداگانه از توزیع‌های حاشیه‌ای 
در نظر گرفت. بنابراین. دشواری یافتن خانواده‌ای از توزیع‌های 
پارامتری که با هر دو توزیع حاشیه‌ای سازگار باشد. مطرح نیست 
(Nelsen, 2006)‏ 

ايدة اصلی تابع مفصل از 4.45 اسکلار است که بیان 
می‌کند برای دو متغیر تصادفی X‏ و ۲ با تابع توزیع توأم 


y)‏ ک ۷ Fyy (x,y) = P(X< x,‏ و توزیع‌های حاشیه‌ای ی 


Fy‏ وجود دارد تابع C‏ بهكونداى که 
Fry (xy) =C(Fr (x). F (»)) 0‏ 


(Sklar, 1959) یک تابع مفصل دومتغیره است‎ € at 
اسکلار‎ dyad عکس‎ 
داریم‎ ) , v) € ]0, 1] x ]0, 1] 


C(u, v) = Fyy (Fx (u), Fy w) و‎ 


که درآ ن © Fx‏ تابع وارون F x‏ و )12 Fy‏ تابع وارون 


در حالت پارامتری C al‏ متعلق به یک خانوادة پارامتری 
از خانوادة مفصل‌های cul ٤= 100:0 EO}‏ که در ol‏ 

al G ab‏ چگالی ,0 است» که در صورت وجود از رابطة 
C (u v) = Cou)‏ به‌دست می‌آید. 


ðuðy 
l برآوردگر حداكثر درست‌نمایی 0 عبارت است‎ 


n 
AML 
ê = arg max ١ logcg (u; vi) 
eO 
i=1 


در جدول ۱ فرمول تابع توزیع و حدود پارامتر چند خانوادة 
معروف از توايع مقصل بیان شذه ایست: 

ضريب همبستگی ۲ کندال با استفاده از تعريف تابع مفصل CC‏ 
به‌صورت زیر است )2006 (Nelsen,‏ 


Txy = Tc = 4 ff C(u, v)dC(u, v) — 1 (0) 


همچنین ضرایب وابستگی دمی بالایی و پایینی» به‌ترتیب. 
عبارت‌اند از )2006 (Nelsen,‏ 


1- Si t) 


2 lim p(y > ۳+ Wx > Fx )0( = =2- in == (۶) 


A = Jim 1 P(Y < ry ox < ۳ )0( = = lim oe Y) 


ضریب وابستگی دمی بالایی (پایینی) بیان می کند که چگونه مقادیر 
بزرگ (کوچک) یک متغیر با مقادیر کوچک (بزرگ) متغیر دیگر رخ 
می‌دهند. خانواده‌های توابع مفصل به ساختارهای وابستگی متفاوتی 
منجر می‌شوند. از جمله اينكه در خانواده‌های پلاکت. نرمال و فرانک 
به‌ازای تمام polis‏ مجاز پارامتره ضرایب وابستگی دمی بالایی و 
پایینی صفر است. در خانواده‌های گامبل. جو و گالامبوس ضریب 
وابستگی دمی بالایی غیرصفر و ضريب وابستگی دمی پایینی صفر 


YVA 


نشربه علمی پژوهشنامه dow‏ دوره AV‏ شماره ۳. تابستان ۱۴۰۳ شماره پیاپی FA‏ ص ۲۸۶-۲۷۵ 


جدول ۱: تابع توزیع و حدود پارامتر خانواده‌های پلاکت. فرانک»گامبل. نرمال. جو و گالامبوس 
Table 1: Distribution functions and the range of their parameter of Plackett, Frank, Gumbel, Normal, Joe and Galambus families‏ 


خانواده حدود پارامتر 
پلاکت }1{-)0,2°( )1 — 6) 4070 — 
فرانک (*,1] 
گامبل ]1,1-[ 
نرمال [1,1-] 
جو («ر1] 


برای متغیرهای تصادفی (T,M)‏ اگر Fr mx (tm |x)‏ 
توزیع توأم آنها به شرط X HX‏ باشد. با توجه به (۳ تابع CL‏ 
نشاندهندة تابع مفصل توأم (T,M)‏ به شرط × = 2 است. در 
اين صورت: 
(tx), Fux nx) (A)‏ رو Fryjx(t, mlx) = Cxe‏ 
در اینجا فرض مىكنيم فقط توزیع‌های حاشیه‌ای (T,M)‏ 
ب‌ستگی ندارد. بنابراین C,‏ = 7 است. 

با فرض V; = Fux M; | X;) 5 U; = Fry (7; |X)‏ تغیرهای 
U;‏ و بلا به‌طور یکنواخت در فاصلة ]1 ,0] توزیع lords‏ و تابع 
توزيع توأم 4 U; ~ Soil G‏ و Ve‏ را به دلیل اينكه 
AN ۳‏ ا وجود 0 cae‏ به‌طور دقيق محاسبه 
کرد. در این صورت پراورد تجربی مقادیر (2 | U, = Fry X;‏ 
و Fy (N, |X;)‏ = محاسبه می‌شوند و برآوردگر پارامتر 
مفصل با توجه به برآوردگر حداكثر درست نمائى وزنى رابطة )\( 
به‌صورت زير خواهد بود: 


۱ 


ô = arg max 2i1 م۱096‎ (û; P;) (a) 
ع6‎ 


W, öl که در‎ 


i‏ در (Y)‏ تعریف شده است. در مورد داده‌های 


سانسورشده. یعنی وقتی 0= ;6 آنگاه 0= ,7 است. بنابراین 


Wa‏ یک وزن برای اصلاح اریبی ناشی از سانسور است. 


پیش‌بینی مقدار نهایی خسارت‌ها 
در قسمت قبل توزیع توأم (T,M)‏ به شرط < 2 با 
استفاده از تابع مفصل ساخته شد. از اين مدل برای پیش‌بینی مقدار 
نهایی هر خسارت استفاذه می‌شود. برای خسارت‌هایی كه سانسور 
نشده‌اند تابع چگالی توام شرطی با استفاده از توزيع توام شرطی 
به‌دست می‌آید. اما برای یک خسارت تسویه‌نشده که برای ol‏ 
0= :۵ است. معلوم است که متغیر مشاهده‌نشدة T,‏ بزرگ‌تر از 
Y,‏ است. بنابراین» توزیع (T,M)‏ به شرط ,× او T2Y,‏ 
مورد نیاز است که دارای برآورد تابع چگالی (tm)‏ بر fy‏ است. 


تابع مفصل 
(0-1)(u+ v) - e + ۷((2‏ +1 = 
a ii D,‏ گ + Co(u,v) = sinc‏ 
Co (u, v)=exp(—((—In u)? + (—1n 8‏ 
(u), P 0((‏ )لوط = Colu v)‏ 
را( 1( °( - 1) - 9( - 1( + Colu v) =1-((1-u)®‏ 


Co(u, v) = uv ١ exp[(—In (+) vy] 


از رابطة (A)‏ تابع چگالی توأم (T,M)‏ به شرط ‏ = 26 عبارت 


fr (tml) = 


(Frx (tx) Fux (m|x)) fry El) Fy mlx)‏ میت 
بنابراین تابع چگالی توأم (TM)‏ به شرط T25X=X‏ 
Y,‏ عبارت است از: 


1 


T,M|X=x,T2Y; 


(£4m) = 
c, o (Fna (t1 x); Eux (m|x)) tnx (tx) fux (m| x)( Ly ) 
اس متأ‎ u | عد‎ Fue (v| |x) Sx (ulx) fine ۳ |x) du dv 


که به‌صورت زیر برآورد می‌شود: 


fx (bm) = 
C; : TIX (t | x), Ê, MIX (m | x)) frx (fx) 7 (m | x) Ley 
= 6 C; (Fax (u u | eae (v | *)) Ax (ule) ane (v | x) du dv 


(+) 


که در آن Cas (BO) x (.|x)‏ تابع چگالی شرطی 
برآوردشده (Fuel |x)) Fa: | x) 9(M|X=x)T|X =x‏ 
تابع توزیع شرطی برآوردشده را نشان می‌دهد. سپس پیش‌بینی مبلغ 
خسارت a)‏ شرط آنکه تا زمان Y,‏ پرداخت نشده باشد) براساس 
براورد اميد ریاضی خسارت از رابطة 


م 00 00+ A‏ 
(t,m)dtdm (\\)‏ الم = E(M|X =x,T >Y,)‏ 
انجام می‌شود. 
شبیه‌سازی 
برای بررسی اعتبار روش پیش‌بینی يادشده در بخش ۲ با مقادیر 


پارامترهای خلاصه‌شده در جدول ۰۲ و سه خانوادة مفصل پلاکت. 


۳۳۹ 


صدیقه شمس و همکاران 


جدول ۲: پارامترهای شبیه‌سازی 
Table 2: Simulation parameters‏ 


نام پارامتر مقدار پارامتر درصد مشاهدات سانسورشده 
01 0 
02 05 
003 02 
0.4 0.1- 
Bo‏ 0.1 
Boa‏ 0.03 
Bo.2‏ 0.05 
Bo.3‏ 0.02- 
(k, A)‏ )2.7,3( 30 
(k, A)‏ )2.5,3( 50 
(k, A)‏ )0.15,3( 60 
جدول ۳: اندازة خطا به‌ازای 0.25 = ۲ و 5000 = n‏ 
Table 3: Error measurement for ۲ = 0.25 and n=5000‏ 
درصد سانسور خانوادة مفصل E, Eir Ey‏ 
30 پلاکت 17.61 0.21 0.014 
فرانک 20.61 0.27 0.012 
28.88 0.30 0.018 
نرمال 
50 پلاکت 15.27 0.15 0.014 
فرانک 12.93 0.27 0.013 
19.41 0.19 0.012 
نرمال 
60 پلاکت 14.04 0.18 0.033 
فرانک 11.04 0.29 0.039 
۱ 17.27 0.26 0.024 
نرمال 


فرانک و نرمال» شبیه‌سازی به شرح زیر انجام می‌شود: 
مرحلة 5 بردارهای 
وک (XX,‏ = ,× شبیه‌سازی می‌شود که در آن bX;‏ 
به‌طور مستقل از توزيع یکنواخت در Abels‏ (۱ و +( توليد می‌شوند. 
مرحلة ۲. به‌ازای 1,...,70 < 7 متغیر T,‏ از توزیع لگ - نرمال 


i=l,...,n به‌ازای‎ 


به صورت شرطی بر روی متغیرهای ,2۶ شبیه‌سازی می‌شود. 
مرحلة ۳. به‌ازای 72و...,1< 7 سانسور D,‏ به‌طور مستقل از 
توزیع وایبل با پارامترهای ۸ و A‏ یعنی با تابع چگالی 
(—d/A)“ Lazo‏ جع ۳0۳ feld) = kA‏ شبیه‌سازی می‌شود. 


٠.5 aby‏ برای هر نمونۀ شبیه‌سازی‌شده 
(X, N, D, 8p X, Jai,‏ برآورد تابع چگالی از رابطة )+( 


مرحلة ۵ به‌ازای 1,...,7= # پیش‌بینی M,‏ از رابطة (VY)‏ 
محاسبه می‌شود. 

Ale yo‏ ۶ مراحل ١‏ تا ۵ را به تعداد 10000 - 8 بار تکرار می كنيم 
تا به‌ازای D=1,...,B‏ نمونه‌های (X, NE, DPO, X?) n,‏ 
شبيهسازى شوند و مقادير براوردشده به دست آیند. 

برای تحليل تأثير قدرت سانسور. Glare‏ مختلفى برای 
سانسور شدن در نظر كرفته می‌شود. اگر =F) (1-6) MP‏ ۸ 
ذخيرة مربوط به تمام مطالبات خسارتها در مرحلة ام باشد و 


e 5 ۰ A ۳ D> _ 4 7‏ 7 3 
4 -1). ,< = ۸ مقدار پیش‌بینی‌شدة ol‏ باشد. 
خطاهای مختلف اندازه‌گیری که به شرح زیر تعریف شده‌اند» محاسبه 


Bia 
1 2 R? — RÈ 

E =— 

17 2 R? 

PE ی‎ 
B i i i 


خطای E,‏ نشاندهندة ميانكين قدر مطلق تفاضل مقدار واقعى 
ذخيره و برآورد آن است. خطای ,رل ميانكين خطای نسبى برآورد 
را اندازه م ىكيرد. خطای Ey‏ ميانكين قدر مطلق تفاضل برآورد هر 
خسارت پرداخت‌نشده و مقدار واقعى خسارت است. 

جدول ۲ نشان مىدهد برای هر Ja‏ خطاى نسبی E,‏ 
تحت توابع مفصل مختلف نسبتا ياياست. وقتى نسبت سانسور 
افزايش می‌یابد pls By‏ به كاهش دارد زيرا تعداد خسارتهاى 


۳۸۰ 


ذخیره خسارت‌های معوق در سطح خرد 


جدول ۴: آمار توصیفی 7 (بر حسب روز) و M‏ (بر حسب ده میلیون ریال) قبل از حذف اثر تورم 
Table 4: Descriptive statistics for T (in terms of day) and M’ (in terms of ten million Rials) before removing the inflation‏ 


متغیر مینیمم چارک اول میانه میانگین چارک سوم ماکزیمم 
M’‏ 0.0027 1.68 5.92 22.14 21.22 715.50 
T‏ 1 85 319 543 828 3234 


پیش‌بینی‌شده افزايش می‌یابد و ازآنجاکه برآوردگرها به‌طور مجانبی 
نااریب هستند. بیش‌برآوردی‌ها را کم‌برآوردی‌ها خنثی می‌کنند. 
خطای By‏ (در مورد برآورد پارامتر مفصل) با نسبت سانسور افزایش 
تایه زرا CI‏ كمقر .يه مدل جا اد قاری متیر 
می‌شود و تأثير منفى بر پیش‌بینی دارد. 


مطالعة موردی دربارة داده‌های واقعى 

پژوهشکدة day‏ یک مجموعه داده با ۲۷۰۲ مورد اعلام‌شده 
ناشی از مسئولیت حرفه‌ای کارفرما در یک دورة ۸ ساله در اختیار 
ما قرار داد که شامل ۲۶۶۰ مورد خسارت بدون سانسور هستند که 
مبلغ نهايى آنها مشخص است 9 ۳۲ مورد سانسور شده است (آنهایی 
که در پایان ذکر است 
كه به‌دلیل نبود داده‌های مناسب که بتوان آنها را به‌عنوان متغیرهای 
کمکی در نظر گرفت. در اين مطالعة موردی از متغیرهای کمکی 
صرف‌نظر شد. هدف برآورد مبلغ خسارت اين FY‏ مورد و برآورد 
ذخیره معوق کل čal‏ 

برای اجتناب از عوامل خارج از کنترل برای تجزیهوتحلیل دادهها 
با روشی که در adol‏ شرح داده می‌شود. عامل تورم حذف مي‌شود. 
متغير 1 (برحسب روز) مدت‌زمان بين اعلام خسا 
است. متغیر M’‏ (برحسب ده میلیون ربال) مبلغ يرداخت شده ازسوى 


اين دورة ۸ ساله تسویه نشده بودند). شایان 


برای موارد تسویه‌شده. نشان می‌دهد . 


حذف اثر تورم 

از روش‌های مختلفی برای اصلاح اثر تورم بر مبلغ خسارت 
استفاده می‌شود. در یکی از ساده‌ترین روش‌ها با توجه به اينكه تورم 
باعث افزایش قيمت‌ها می‌شود. فرض می‌شود افزايش ساليانة قيمت 
ثابت باشد و لگاریتم نسبت قيمت‌ها در دو سال متوالی ضریبی از 
سال فرض شده است. به‌عبارت‌دیگر: 


M, 
= HS: 
M, 
مقدار‎ M, مقدار خسارت را قبل از حذف تورم‎ M, كه در آن‎ 
در‎ ali خسارت بعد از حذف تورم. ,© سالی است که خسارت‎ 
€ 10,1,2,..., G} آن اتفاق افتاده است و فرض مىكنيم‎ 
ضريبى ثابت است که با استفاده از مجموعه داده‌ها‎ 1 
مستقل از‎ (MG) icy برآورد می‌شود. در اين مدل فرض می شود‎ 
بر (:8) است. و بنابراین از مدل زیر استفاده می‌شود.‎ 


و ضریب 


log M; = ۸.9: + log M; ۳‏ 
برای برآورد دا روش زیر به کار می‌رود: 

* قرار می‌دهیم M;‏ را به‌عنوان میانگین تعداد Ni‏ خسارت‌هایی 
که در سال Jh‏ ,£ رخ داده است و يس از 7 سال تسویه شده 
است؛ 

«. از رابطة +A;‏ ,لم + 1087 که در ol‏ 
Ay = Eflog M; |T; = j]‏ برای هر J‏ است. برآوردگر وزنی 

حداقل توان دوم خطا را برای ,لم محاسبه م ىكنيمء یعنی: 


2 
arg mın dig, +i<G ni (log mj — a — bg;) 


* برآوردگر نهایی ۸۸ عبارت است از: 


(VY) 


كه در آن 20 = ر. 

در نهایت» برای هر i‏ كه 1= :8 است مقدار خسارت پس از 
حذف اثر تورم به صورت log M, = log M; - fig,‏ محاسبه مىشود. 
جدول ۵ آمار توصيفى متغير M‏ بعد از حذف اثر تورم را نشان 
می‌دهد. 

نمودارهای ۱ و ۲ هیستوگرام مربوط به مبالغ خسارات 
پرداخت‌شده (بعد از حذف اثر تورم) و هیستوگرام مدت‌زمان پرداخت 
خسارات را نشان می‌دهند. 


برآورد توزیع‌های حاشیه‌ای 
با توجه به نمودارهای هیستوگرام متغیرهای 7 و ۸ که 
حا از جولكى ميت ای كالى: هر کو تفر veal‏ تور ای 
نمایی» گاماء لگ نرمال و dels‏ را برای برازش به آنها آزمون کردیم؛ 
ولی هیچ‌یک از آنها برازش مناسبی نداشتنده لذا تصمیم گرفتیم از 
تقریب توزیع نرمال استفاده کنیم. مناسب‌ترین تبدیل برای متغیر 
M‏ تبدیل, و Box-Cox‏ با ضریب 0.0554 = A‏ به‌صورت 
M“ =h (m) - ,A#0‏ و مناسب‌ترین تبدیل 
متغير pos T‏ معروف Order- Norm‏ به‌صورت h, (t)‏ = 
(04-(۳۷ $ = تشخیص داده شد. جدول ۶ نتایج os‏ به 
آزمون نرمالیتی را برای متفیرهای "7 و M`‏ نشان مىدهد. 
نمودارهاى 7 و ۴ هيستوكرام مربوط T”‏ و M?‏ هستند. 


YA\ 


نشریه ole‏ پژوهشنامه dow‏ دوره AT‏ شماره ۰۳ تابستان ۰.۱۴۰۳ شماره پیاپی ۰۴۹٩‏ ص ۲۸۶-۲۷۵ 


جدول ۵: آمار توصیفی M‏ (برحسب ده میلیون ریال) بعد از حذف اثر تورم 
Table 5: Descriptive statistics and M (in terms of ten million Rials) after removing the inflation‏ 


متغیر مینیمم چارک اول میانه 
M‏ 0.0022 1.15 4.07 
an‏ كبر a‏ 


میانگین چارک سوم ماکزیمم 
16.03 15.20 517.33 
ل | ل هه 
Š 8 pa‏ 


Claim amount (10 Million Rials) 


نمودار ۱: هیستوگرام مبلغ خسارت. 
Fig. 1: Histogram of the claim amount‏ 
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نمودار ۲: هیستوگرام مدت‌زمان پرداخت مبلغ خسارت. 
Fig. 2: Histogram of the duration of the payment of claim‏ 


* 
جدول ۶: نتایج مربوط به آزمون نرمالیتی را برای متغیرهای ‏ و 


* 


M 


Table 6: Results of the Normality test for T and M* 


متغیر آمارة آزمو ن کلم وگ روف اسمیرنوف مقدار احتمال 
M*‏ 0.0168 0.4382 
T*‏ 0.0101 0.9469 


برای بیان وابستگی مبلغ خسارت پرداخت‌شده (پس از تبدیل به 
نرمال) در مقابل مدت‌زمان پرداخت خسارت (پس از kas‏ به نرمال) 
نمودار پراکنش آنها رسم شده است نمودار ۵ و همچنین برآورد 
نمونه‌ای ضريب وابستكى 1 كندال و ضرايب وابستكى دمى بالايى و 
پایینی را برای داده‌های متغيرهاى T‏ و M‏ محاسبه شده‌اند که 


در جدول ۷ نشان داده شده است. 


برازش تابع مفصل ۱ 

all‏ مفصل تحت تبدیل‌های اکیدا افزایشی. پایاست بنابراین 
تابع مفصل متفیرهای "7 و 17 همان تابع مفصل متغيرهاى 
(7,M)‏ است. از بين توابع مفصل مختلفی که به داده‌های زوج 
M)‏ , 7) برازش داده شد. براساس آزمون نیکویی برازش, فرضية 
صفر برازش fare al jhe‏ پلاکت. فرانک» نرمال و جو در سطح 


YAY 


صديقه شمس و همکاران 


1 1 1 1 
0 ۹ 2 3 


Frequency 


r T T 
-3 2 27 


Normal Transformed of Claim amount 


نمودار iY‏ هیستوگرام مبلغ خسارت پس از تبدیل به نرمال. 
Fig. 3: Histogram of the claim amount after Normal transformation‏ 
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نمودار ۴: هیستوگرام مدت‌زمان پرداخت ale‏ خسارت يس از تبدیل به نرمال. 
Fig. 4: Histogram of the duration of the payment of claim after Normal transformation‏ 


معناداری 0.01 رد نشد. جدول A‏ مقدار احتمال آزمون نیکویی 
برازشء مقدار لگاریتم تابع درست‌نمایی؛ برآورد پارامتر و برآورد 
ضریب ۲ کندال برای هریک از اين توابع مفصل را نشان می‌دهد. 

برآورد نمونه‌ای ضرایب وابستگی دمی براساس داده‌ها با هم 
برابر و تقريباً برابر صفر است. بنابراین خانواده‌های پلاکت. نرمال و 
فرانک گزینه‌های مناسبی هستند و تابع مفصل جو از اين نظر مناسب 
نیست: از بین اين سه ختانواذه.براسان ماک مقدار انم :درست تمان ge‏ 
چون لكاريتم تابع درست‌نمایی تابع مفصل بلاكت از بقيه بيشتر است 
تابع مفصل يلاكت را به‌عنوان كزينة نهايى تابع مفصل برای زوج 
( 26 , 7) انتخاب مى كنيم. 


برآورد ذخيرة خسارت معوق 

با برازش تابع مفصل پلاکت به ( (TM‏ با استفاده از رابطة 
ab )۱۰(‏ چگالی توأم بددست می‌آید و برآورد خسارت مواردی که 
تسویه نشده‌اند با استفاده از رابطة (۱۱) انجام می‌شود. تبدیل وارون 
باکس کاکس روی M 5 polis‏ اعمال شده و در Cols‏ برآوردی 
از مقدار ذخیره به‌صورت ,)1-8( R=Y‏ به دست می‌آید. در 
داده‌های مورد مطالعه ۴۲ مورد تسويەنشدە در بين ۲۷۰۲ خسارت 


اعلام‌شده وجود داشت که برآورد شده و نتایج در جدول ٩‏ نشان داده 
شده است. 

برای ارزیابی خطای برآورد با اين روش در اين مجموعه داده. 
از بين داده‌هایی که ale‏ خسارت انها معلوم است و تسویه شده به 
تعداد ۱۰۰۰ بار نمونه‌هایی به حجم ۴۰ (هر بار فرض شده که ۴۰ 
مورد تسویه‌نشده وجود دارد) شبیه‌سازی کردیم. در اين نمونه‌ها مبلغ 
خسارت معلوم است. ولی اين مبالغ برآورد نيز شده‌اند و خطاهای 
E,,E,, , E,‏ 
است. 

نتایج جدول ۱۰ حاکی از آن است که در 
نمونه‌های ۴۰تایی از خسارت‌های eagles‏ میانگین مقدار واقعی کل 
خسارت پرداخت‌شده ۴۰ مورد برابر ۳۲۱/۲۱ (ده میلیون ریال) است 
و با فرض اينكه برای اين نمونه‌ها مبلغ خسارت سانسور شده است. 
میانگین برآورد کل خسارت برابر ۴۹۸,۹۱ «ده میلیون (JL‏ متوسط 
قدر مطلق خطای برآورد خسارت کل برابر ۲۶/۳۹ (ده میلیون ریال). 
متوسط خطای نسبی برآورد خسارت کل برابر ۰/۱۲۹ و ميانكين 
میانگین قدر مطلق هر مورد خطا برابر ۰/۷۹ (ده میلیون ریال) است. 
در نمودار ۶ هیستوگرام فراوانی برآورد ذخيرة کل به‌دست‌آمده در 


م‌حاسبه و نتایج در جدول ٩‏ نشان داده شده 
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نمودار ۵: نمودار پراکنش مبلغ خسارت در مقابل مدت‌زمان پرداخت خسارت يس از تبدیل به نرمال. 
Fig. 5: Scatter plot of claim amount versus the duration of its payment after Normal transformation‏ 


جدول ۷: برآورد نمونه‌ای از ضریب وابستگی ۲ کندال و ضرایب وابستگی دمی بالایی و پایینی 
Table 7: Sample estimations of the Kendall’s ۲ and upper and lower tail dependence coefficients‏ 


Àu 2 T زوج متغیر‎ 
005 005 0.25 T*,M* 


جدول ۸: نتایج آزمون نیکوبی برازش و برآورد پارامترهای توابع مفصل منتخب 
Table 8: Results of the good ness of fit test and the parameter estimations of the selected copula function‏ 


خانوادة مفصل مقدار احتمال لکاریتم درست‌نمایی بر آورد پارامتر برآورد 7 بر آورد (Ar, Ay)‏ 
پلاکت 0.0528 5.92 0.0762 0.2455 (0,0) 
نرمال 0.0544 5.45 0.3455 0.2232 )0,0( 
فرانک 0.0244 5.23 2.1043 0.2241 (0,0) 
جو 0.1043 7.82 1.4791 0.2122 (0,0.04) 


جدول :٩‏ برآورد مبلغ خسارت‌های تسویه‌شده (برحسب ده میلیون ریال) 
Table 9: Estimations of unsettled claim amount (in terms of ten million Rials)‏ 


ردیف برآورد مبلغ ردیف بر آورد مبلغ Wo,‏ برآورد مبلغ ردیف برآورد مبلغ 
10 8.63 772 14.79 1766 10.65 2680 9.04 
166 10.57 940 9.28 2660 8.87 2683 7.99 
167 10.64 941 13.76 2665 10.53 2684 9.13 
345 17.59 963 13.40 2666 8.77 2685 8.57 
445 12.12 1453 12.51 2668 11.28 2687 9.46 
446 12.20 1454 12.13 2671 9.71 2692 10.03 
569 14.16 1455 9.67 2672 10.28 2696 9.49 
648 11.38 1456 12.85 2674 9.31 2698 9.93 
649 11.47 1457 13.45 2675 9.51 2700 9.94 
770 14.37 1514 9.67 2676 10.01 جمع 
R= 463.50 10.74 2678 11.03 1595 14.50 771‏ 


جدول ۱۰: متوسط مقدار واقعى خسارات. متوسط مقدار برآوردشدة ذخيره و خطاهاى برآورد در ٠٠٠١‏ بار شبیه‌سازی 
Table 10: The average of real amount of claims, the average of estimated amount of claims and estimation errors in 1000 times‏ 
of simulation‏ 


Ez و‎ A R R 
0.79 0.129 26.39 498.91 481.21 
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نمودار ۶: فراوانی برآورد ذخیره خط قرمز نشان‌دهندة میانگین واقعی ذخیره 
Fig. 6: Reserve estimate frequency, red line represents the actual average reserve‏ 


شبیه‌سازی و میانگین مقدار واقعی ذخیره (خط قرمزرنگ) olis‏ 


داده شده است. 


جمع بندی و پيشنهادها 

در اين ale‏ برای پیش‌بینی مبلغ نهایی یک خسارت 
تسویه‌نشده» از ساختار وابستگی بين مقدار خسارت و طول عمر ol‏ 
استفاده شد. اين ساختار وابستگی با تابع مفصل که امکان توصیف 
وابستگی را به‌صورت جداگانه‌ای از مدل توزیع‌های حاشیه‌ای فراهم 

برای بهتر شدن دقت پیش‌بینی در صورت موجود بودن 
خصوصیات شناخته‌شده‌ای از مبلغ خسارت و طول عمر آن (به‌عنوان 
متغیرهای کمکی) می‌توان مدل‌سازی حاشیه‌ای را براساس توزیع 
شرطی به اين خصوصیات انجام داد. برای انتخاب تابع مفصل مناسب 
به زوج متغیر مبلغ خسارت و طول عمرء از تطابق برآورد ضرایب 
وابستگی دمی نمونه‌ای ( ناپارامتری) و برآورد ضرایب وابستگی دمی 
ali‏ مفصل استفاده کردیم. از بين توابع مفصل‌های انتخاب‌شده در 
مرحلة تطابق, با استفاده از آزمون نیکویی برازش توابعی که در سطح 
معناداری ١‏ ۰,:رد نمی‌شدند گزینش شدند و سپس با معیار بیشترین 
مقدار ali‏ درست‌نماییء al‏ مفصل پلاکت که بیشترین مقدار 
درست‌نمایی را داشت به‌عنوان ali‏ مفصل منتخب برگزیده شد. يس 
از انتخاب تابع مفصل با استفاده از قضية اسکلار تابع توزیع و تابع 
چگالی دومتغیره مبلغ خسارت و مدت‌زمان پرداخت به دست آمد. 
ازآنجاکه برای هر خسارت تسویه‌نشده» مدت‌زمان پرداخت از زمان 
شروع مطالعه بيشتر است تابع چگالی شرطی مبلغ آن خسارت به 
شرط آنکه مدت‌زمان پرداخت ان از يك مدت مشخص بیشتر باشد. 
محاسبه شد. سپس مبلغ هر خسارت تسویه‌نشده با استفاده از تقریب 
میانگین شرطی أن پیش‌بینی شد و جمع مقادیر اين پیش‌بینی‌ها 
به‌عنوان پیش‌بینی کل خسارت معوق در نظر گرفته شد. همان گونه 
که شرح داده شد در اين روش انتخاب تابع مفصل و انتخاب متغیرهای 
کمکی اهمیت بسیار زیادی در دقت پیش‌بینی دارند. 


شانان Sail oy‏ که انتفامه او gal‏ رون لمات خا ا 
اطلاعات مورد نیاز با دقت ثبت و نگهداری شوند. 


صدیقه شمس: طرح موضوع. ارائة «fool,‏ تدوین alio‏ 
ارائة الگوریتم محاسباتی» تنظیم کدهای برنامة کامپیوتری» مریم 
اننی‌عشری: پردازش مسئله. مشاوره در بسط «fool,‏ يافتن داده‌هاء 
کمک در تدوین alio‏ محبوبه پیاده کوهسار: پالایش ۳ تنظیم 
داده‌هاء کمک در نوشتن کدهای برنامه کامپیوتری. 


از پیشنهادهای داوران محترم که به غنای علمی مقاله کمک 
کردند. بسیار سپاسگزاريم. 
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